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Abstract 



We present a solvable two-dimensional piecewise linear chaotic map which arises from the 
duplication map of a certain tropical cubic curve. Its general solution is constructed by means 
of the ultradiscrete theta function. We show that the map is derived by the ultradiscretization 
of the duplication map associated with the Hesse cubic curve. We also show that it is possible 
to obtain the nontrivial ultradiscrete limit of the solution in spite of a problem known as "the 



minus-sign problem." 

1 Introduction 

> . 

Ultradiscretization [41] has been widely recognized as a powerful tool to extend the theory of in- 
tegrable systems to the piecewise linear discrete dynamical systems (ultradiscrete systems) [7, 11, 
23,29, 31,36, 37,39,44,47,50]. In particular, it yields various soliton cellular automata when it is 
possible to discretize the dependent variables into finite numbers of integers by a suitable choice of 
parameters. It has also established the links between the theory of integrable systems and various 
areas of mathematical sciences, such as combinatorics, representation theory, tropical geometry, 
traffic flow models, and so on [5,6, 10, 13, 17, 19-22,27,30,31,33,43,45,46,49]. One of the key 
features of ultradiscretization is that one can obtain piecewise linear discrete dynamical systems 
from rational discrete dynamical systems by a certain limiting procedure, which corresponds to 
the low-temperature limit in the statistical mechanics. When this procedure is applied to a certain 
class of discrete integrable systems, wide classes of exact solutions, such as soliton solutions or pe- 
riodic solutions, survive under the limit, which yield exact solutions to the ultradiscrete integrable 
systems. 

While the application of ultradiscretization to the integrable systems has achieved a great suc- 
cess, it seems that only few results have been reported as to application to non-integrable sys- 
tems [38,40]. One reason may be that in many cases of non-integrable systems, fundamental 
properties are lost under the limit. For example, it is possible to ultradiscretize the celebrated 
logistic map formally, but however, its chaotic behavior is lost through the ultradiscretization [16]. 
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In [16], the ultradiscretization of a one-dimensional chaotic map which arises as the duplication 
formula of Jacobi's sn function has been considered. It exemplifies a solvable chaotic system, 
which is regarded as a dynamical system lying on the border of integrability and chaos [4], in the 
sense that though its exact solution is given by an elliptic function, however, its dynamics exhibits 
typical chaotic behaviors such as irreversibility, sensitivity to the initial values, positive entropy, 
and so on. By applying the ultradiscretization, it has been shown that we obtain the tent map and 
its general solution simultaneously. Moreover, a tropical geometric interpretation of the tent map 
has been presented, namely, it arises as the duplication map on a certain tropical biquadratic curve. 
This result implies that there is the world of elliptic curves and elliptic functions behind the tent 
map, which might be an unexpected and interesting viewpoint. It also suggests that the tropical 
geometry and the ultradiscretization provides a theoretical framework for the description of such a 
geometric aspect. 

In this paper, we present two kinds of two-dimensional solvable chaotic maps and their general 
solutions that are directly connected through the ultradiscretization. In Section 2, we construct a 
piecewise linear map from a duplication map on a certain tropical plane cubic curve. We also con- 
struct its general solution in terms of the ultradiscrete theta function [16,20,26,32,33,42] by using 
the tropical Abel-Jacobi map. In Section 3, we consider a certain rational map which arises as a 
duplication map on the Hesse cubic curve (see, for example, [1, 15,35]), whose general solution 
is expressible in terms of the theta functions of level 3. In Section 4, we discuss the ultradis- 
cretization of the rational map and its solution obtained in Section 3, and show that they yield the 
piecewise linear map and its solution obtained in Section 2. The rational map and its general solu- 
tions discussed in Section 3 and 4 involve a problem known as "the minus-sign problem," which is 
usually regarded as an obstacle to successful application of the ultradiscretization. We show that it 
is possible to overcome the problem by taking careful parametrization and limiting procedure. 

2 Duplication map on tropical cubic curve 

2.1 Duplication map 

In this section, we construct the duplication map on a certain tropical curve. For basic notions of 
the tropical geometry, we refer to [3, 12, 24-26, 34]. 

Let us consider the tropical curve C K given by the tropical polynomial 

x ¥(X,Y;K) = maz[3X,3Y,X + Y + K,Q], X,Y,KeR, K>0. (2.1) 

The curve C K is defined as the set of points where *P is not differentiable. As shown in FigOJa), 
the vertices V; and the edges E t of C K are given by V\ = (-K, 0), V 2 = (0, -K), V 3 = (K, K) and 
Ei = V1V2, E 2 = V2V3, £3 = V3V1, respectively. From the Newton subdivision of the support of 
Cr given in FigOJb), we see that Ck is a degree 3 curve. For a vertex on the tropical curve, let 
v, e Z 2 (z = 1, . . . ,ri) be the primitive tangent vectors along the edges emanating from the vertex. 
Then it is known that for any vertex there exist natural numbers w, e Z >0 (i = 1, . . . ,n) such that 
the following balancing condition holds: 

wivi + --- + w„v n = (0,0). (2.2) 
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We call Wj the weight of corresponding edge. Now, since the primitive tangent vectors emanating 
from V\ is given by (-1,0), (1, -1) and (2, 1), the balancing condition at V\ is given by 



3(-l,0) + (1,-1) + (2,1) = (0,0). 



(2.3) 



Therefore, the weight of E\, Zs 3 and the tentacle along the edges emanating from V\ are given 
by 1, 1 and 3, respectively. The balancing condition at V 2 and V 3 shows that the weight of the 
edges Ej (i = 1,2,3) are 1, and those of tentacles of C K are all 3, respectively. If a vertex V is 
3-valent, namely V has exactly three adjacent edges whose primitive tangent vectors and weights 
are v, and w t (i = 1,2,3), respectively, the multiplicity of V is defined by w\w 2 |det(vi, v 2 )| = 
W2W3 |det(v2, v 3 )| = W3W1 |det(v 3 , Vi)|. If all the vertices of the tropical curve are 3-valent and have 
multiplicity 1, then the curve is said to be smooth. The multiplicity of the vertex V\ is computed 
as 

-1 
1 -1 



3-1 



det 



= 3, 



(2.4) 



and similarly those of V 2 and V 3 are both 3, which imply that C K is not smooth. The genus is equal 
to the first Betti number of C K , which is 1 as shown in FigfjTa). Thus the curve C K is a non-smooth, 
degree 3 tropical curve of genus 1 . Note that the cycle Ck of Ck (the triangle obtained by removing 
the tentacles from Ck) can be given by the equation 



max[3X,3y,0] = X + Y + K. 



(2.5) 




(a) (b) (c) 

Figure 1. (a): Tropical curve C K . V\ = (-^,0), V 2 = (0,-K), V 3 = (K,K). Primitive tangent 
vector of each edge: E\. v x = (1,-1), E 2 : v 2 = (1,2), E 3 : v 3 = (2, 1). (b): Newton subdivision of 
the support of C K . (c): Tropical line. 



A tropical line is the tropical curve given by the tropical polynomial of the form 

L(X, Y) = max[X +A,Y + 5,0], (2.6) 

which is shown in Fig{T](c). The three primitive tangent vectors emanating from the vertex are 
given by (-1, 0), (0, -1), (1, 1). From the balancing condition (-1, 0) + (0, -1) + (1, 1) = (0, 0), the 
weight of edges are all 1 . 

Vigeland [48] has introduced the group law on the tropical elliptic curve, which is a smooth, 
degree 3 curve of genus 1. According to the group law, the duplication map is formulated as 



3 



follows; let C be a tropical elliptic curve and let C be its cycle. Take a point P e C. We draw a 
tropical line that intersects with C at P with the intersection multiplicity 2, and denote the other 
intersection point by P * P. Drawing a tropical line passing through O and P * P with a suitable 
choice of the origin of addition O e C, the third intersection point is IP. 

For a given point P on a tropical elliptic curve, the tropical line that intersects at P with the 
intersection multiplicity 2 does not exist in general. However, the curve C K has a remarkable prop- 
erty that it is possible to draw a tropical line that intersects at any point on C K with the intersection 
multiplicity 2. The explicit form of the duplication map is given as follows: 

Proposition 2.1 Choosing the origin as O = V3, the duplication map Ck 3 P = (X, Y) 1 — > 2P = 

(X, Y) 6 C K on the tropical cubic curve C K is given by 

X = Y + 3 max[0, X] - 3 max[X, Y], Y = X + 3 max[0, Y] - 3 max[X, Y], (2.7) 



or 



X n+1 = Y n + 3 max[0, X n ] - 3 max[Z„, Y n ], Y n+1 =X n + 3 max[0, Y n ] - 3 max[X„, Y n ], (2.8) 
where (X n , Y n ) is the point obtained by the n times successive applications of the map to (X, Y). 
Proof. 

Case 1: P € Ei As illustrated in Figf2|a), the primitive tangent vectors of the two edges 
passing through P are (1,-1) (thick line) and (1,1) (broken line), respectively, and the weight of 
the edges crossing at P are both 1. Then the intersection multiplicity is given by 



1 • 1 • 



det 



1 -1 
1 1 



2. 



(2.9) 



Let P * P be the other intersection point. Note that the intersection multiplicity at P * P is 1 . Then 
the map P = (X, Y) 1 — > P*P = (X', Y') is constructed as follows. Since P e E Y and P*P e Ei UE 2 , 
we have 

= X + Y + K, max[3X', 3F] = X' + Y' + K. (2.10) 
Subtracting the second equation from the first one, we obtain by using (X - X')/(Y -Y') = \ 

X' = X - 3 max[X, Y], T = Y - 3 max[X, Y], (2.1 1) 

Our choice of the origin of addition O = V3 makes the form of IP simple. It is obvious as illustrated 
in Figf2](b), that IP = (X, Y) is given by (X, Y) = (Y',X'). Hence we obtain the map P 1— > IP as 



X = Y - 3 max[X, Y], Y = X - 3 max[X, Y], (X, Y) 6 E x . 



(2.12) 



Case 2: P € E 2 As illustrated in FigfJta), the two primitive tangent vectors of the edges 
passing through P are (1,2) and (1,0), and the weight of the edges crossing at P are both 1. Then 
the intersection multiplicity is given by 



1 • 1 • 



det 



1 2 
1 



(2.13) 
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Figure 2 . (a) : Map P i-» P*P for P & E\. The intersection point of CV and the line pas sing through 
P with multiplicity 2 (broken line) is P * P. (b): Map P * P i-» 2P. The intersection point of C# 
and the line passing through (9 = V 3 and P * P (broken line) is IP. Obviously P * P and IP are 
symmetric with respect to X = Y. 




(a) (b) 
Figure 3. (a): Map P i-» P * P for P e E 2 . (b): Map P*P^2P. 



Since P e E 2 and P * P e Ei U £3, we have 

3X = X + y + Z, max[3F\ 0] = X' + Y' + K. (2.14) 
By using Y' = Y, we obtain the map P 1 — > P * P and P 1 — > 2P as 

X' = -2X + 3 max[0, Y], T = Y, (2.15) 

X = Y, Y = -2X + 3 max[0, Y], (X, Y) e E 2 , (2.16) 

respectively. 

Case 3: P € £3 As illustrated in FigHJa), the two primitive tangent vectors of the edges 
passing through P are (2, 1) and (0, 1), and the weight of the edges crossing at P are both 1. The 
intersection multiplicity is given by 



1 • 1 



det 



2 1 
1 



2. (2.17) 



Since P e E 3 and P * P e E\ U E 2 , we have 

3Y = X + Y + K, max[3X', 0] = X' + Y' + K. (2.18) 
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v 2 \ 

(a) (b) 
Figure 4. (a): Map P i-> P * P for P e E 3 . 



(b): Map P * P i-» 2P. 



By using X' = X, we obtain the map P i-» P * P and P i-> 2P as 

X' = X, F' = -2Y + 3max[0,X], (2.19) 
X = -2Y + 3 max[0, X], Y = X, (X, Y) e E 3 , (2.20) 

respectively. 

We finally obtain eq. (12771 ) by collecting eqs. (12.121 ), (12.161) and (12.201 ) together. □ 

Strictly speaking, the group law in [48] cannot be applied to our case, since C K is not a smooth 
curve. However, it is possible to show by direct computation that the map (|2.8I) is actually a 
duplication map on the tropical Jacobian of Ck- For this purpose, we first compute the total lattice 
length X of Ck, which is defined by the sum of the length of each edge scaled by the norm of 
corresponding primitive tangent vector: 

>n \Ei\ V5X V2X V5X 
£ = V^ = _ + -I_ + JL- = 3X. (2.21) 

ti ' v 'l V5 V2 V5 
Then the tropical Jacobian J(C K ) of C K is given by 

J(C K ) = R/XZ = R/3KZ. (2.22) 
The Abel-Jacobi map ?] : C K — » J(C K ) is defined as the piecewise linear map satisfying 

rfiO) = r](V 3 ) = 0, rj(V,) = l ^4=K, rj(V 2 ) = i/(Vi) + ^ = 2K. (2.23) 

|v 3 | |v 2 | 



Proposition 2.2 T/ze map C K 3 P = (X, Y) i — > P = (X, Y) e C K defined by eq. A2.7\) is a 
duplication map on the Jacobian J(Ck)- Namely, we have r](P) = 2r](P) mod 3K. 

Proof. We consider the case P 6 E\. Suppose P = (X, Y) satisfies V\P : V 2 P = s : 1 - s 
(0 < s < 1), namely 

P = (X, Y) = (-(1 - s)K, -sK), 7/(P) = rjiVO + sK = (l + s)K. (2.24) 
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Case (I): X < Y (0 < s < ±). From eq. (12TT21) . P is given by 



X = Y-3Y = -2Y = 2sK, Y = X - 3Y = (-1 + 4s)K, P e E 2 , (2.25) 

which implies 

V 2 P : V 3 P = 2s : 1 - 2s, 77(F) = t](V 2 ) + 2sK = 2(1 + s)K = 2?](P). (2.26) 
Case (II): X > Y (j < s < 1). In this case, P is given by 

X=Y-3X = (3-4s)K, Y = -2X = 2(1 - s)K, P e E 3 , (2.27) 

which implies 

V 3 P : V l P = -1 +2s : 2(1 - s), 77(F) = (-1 + 2s)K = 2(1 + s)K = 2r](P) mod 3K. (2.28) 

Therefore we have shown that tj(P) = 2rj(P) for P e Ei. We omit the proof of other cases since 
they can be shown in a similar manner. □ 

2.2 General solution 

From the construction of the map (12.81) . it is possible to obtain the general solution by using the 
Abel-Jacobi map of C K . Let n\ and n 2 be projections from C K to the X-axis and the F-axis, 
respectively. Then the maps n\ o rf x and n 2 o if 1 , namely, the maps from the tropical Jacobian to 
the X-axis and the Y-axis through the Abel-Jacobi map are given as illustrated in Fig. [5£a) and (b), 
respectively. Therefore, Proposition 12.21 implies that X n = n\ o T]~ l (2 n u ), Y n = n 2 o i]~ 1 (2 n u ) for 



X Y 




arbitrary w 6 J(C K ) gives the general solution to eq. (12.81) . 

It is possible to express n\ o rf l and n 2 o vf x by using the ultradiscrete theta function @(u; 6) 
defined by [16,20,26,32,33,42] 

0(w; 6) = -6 j((w)) - i j , ((«)) = u - Floor («). (2.29) 
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00* -0), 



S{u;a,0, 9) 











0-a \ 


/ P 







00* -0), 



Figure 6. Graph of S{u; a, ft, 6). 



For this purpose, we introduce a piecewise linear periodic function S («; a,fi, 6) by 

<u-fi 



S (u; a, (5, G) = fl) - 0) , 



(2.30) 



which has a period or and amplitude 2/3(a - /3)6/a 2 as illustrated in Figj6l Comparing Fig. [5] with 
Fig. [6l we have 

«^ = s(u-K„\k), W(«) = 4-2™,?4 (2.31) 
Therefore, we obtain the following proposition: 



Proposition 2.3 For a given initial value Pq = (Xq, Yq), the general solution to the map ( 12.ffD is 

given by 

X n = s{l n u -K;3K,2K 9 -K^, Y n = S {tuq - IK; 3K, K, ^Kj , 
K = 3 max[X , Y , 0]-X - Y , u = n(P ). 



Fig. |7] shows the orbit of the map (12.81 ) plotted with 3,000 times iterations. The map has an 
invariant curve C K given by eq. (12.51) . and the figure shows that the curve is filled with the points 
of the orbit. 



3 Duplication map on Hesse cubic curve 
3.1 Duplication map 

The Hesse cubic curve is a curve in P 2 given by 

: x 3 + y 3 + 1 = 3/j.xy, 
or in the homogeneous coordinates [x : X\ : x 2 ] = [x : y : 1] 

Ejj '. Xq + x 3 + x 3 ) = 3\xxqX\X 2 . 



(3.1) 



(3.2) 
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The nine inflection points are given by [1 : -1 : 0], [1 : -co : 0], [1 : -co 2 : 0], [1 : : -1], 
[1:0: -to], [1:0: -co 2 ], [0:1: -1], [0:1: -co], [0:1: -co 2 ], where co is a nontrivial third 
root of 1. It is known that any non-singular plane cubic curve is projectively equivalent to E M (see, 
e.g. [1]). Moreover, these inflection points are also the base points of the pencil 

t Q (x 3 + x\ + x\) = hx xix 2 , [t : t { ] e P 1 . (3.3) 

The duplication map is constructed by the standard procedure; for an arbitrary point on P e E M 
draw a tangent line, and set the other intersection of the tangent line and E M as P * P. Taking one of 
the inflection points as an origin O of addition, the intersection of E M and the line connecting P * P 
and O gives IP. Choosing O to be [1 : -1 : 0] among nine inflection points of E M , the duplication 
map P = (x,y) i-» 2P = (x,~y) is explicitly calculated as (see, for example, [15,35]) 

_ = o-^ (1^ 

x 5 - y i y i — x s 

or writing the point obtained by the n times applications of the map to (x,y) as (x n ,y n ), we have 

_ _ (i -y^xn 

%n+l — 3 3 j yn+l — 3 3 ■ W-Jj 
X n ~ yn yn ~ X n 

By construction, it is obvious that the map (13.51) has the invariant curve E M , where /j. is the conserved 
quantity. Fig. [F] shows the orbit of the map (13.51) plotted with 3,000 times iterations. Note that, 
although the invariant curve has a component in the first quadrant x,y > for jx > 0, the real 
orbit never enters in this quadrant (except for the initial point), which can be verified by a simple 
consideration; suppose that x n > at some n. Then eq. (13.51) implies that (x n _i,j n _i) must be in 
the highlighted region of Fig. H^a). On the other hand, if y n > at some n, (x n -\,y n ~\) must be 
in the highlighted region of Figure 1Mb). Since the intersection of the two regions is empty, it is 
impossible to realize x n ,y n > for any n as long as we start from the real initial value. 
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Figure 8. Orbit of the map (13.51 ) with the initial value (x ,y ) = (2.1,5.3). Dashed line of the left 
figure is the invariant curve. 



y ' 








= X 

X 




l 





(a) (b) 

Figure 9. (a): Region of (x n ^,y n -i) for x„ > 0. (b): Region of (x n - U y n -i) for y n > 0. 

3.2 General solution 

The general solution to the map (13.41) or (13.51) is given in terms of the following theta functions of 
level 3. Let us introduce the functions 6t(z, r) (k = 0, 1, 2) by 

B k { z , r) = Y i e'^-^e 6 ^-^^ = %-i,|)(3z, 3r), (3.6) 

where #( Ui /,)(z, r) is the theta function with characteristic (a, b) defined by 

${a,b)(z, t) = Yj e ni{n+a)lT+2ni{n+a){z+h \ t e <H = {Im z > 0, z e Q. (3.7) 

Proposition 3.1 The general solution to eq. ( 13.51) is given by 

_ e (Tz ,T) _ e 1 (2"z ,T) 

6 2 (2 n zo, r) 2 (2 n zo, r) 

where zo 6 C w an arbitrary constant. 

Proposition [3J] is a direct consequence of the following proposition: 
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Proposition 3.2 

(1) 9 k (z,r) (k = 0,1, 2) satisfy 

9 (z, r) 3 + 9i(z, r) 3 + 2 (z, r) 3 = 3//(t) o (z, t)0i(z, t)9 2 (z, t), /u(t) 
where 

(2) 9u{z,t) (k = 0, 1,2) satisfy the following duplication formulas: 

9 (0, t) 3 9 (2z, r) = r) [fl 2 (z, r) 3 - o (z, r) 3 ] , 
O (O, r) 3 1 (2z, r) = o fc r) [0i (z, r) 3 - 2 (z, r) 3 ] , 
O (O, r) 3 2 (2z, r) = 9 2 (z, r) [0 o (z, r) 3 - 9 t (z, r) 3 ] . 

It seems that the above formulas are well known [2], but however, it might be useful for non- 
experts to give an elementary proof here. In the following, we fix r e 'H and write 9^(z, t) = 9k(z). 

Lemma 3.3 9k(z, t) (k = 0, 1, 2) satisfy the following addition formulas: 



9 (.0) 2 9 (x + y)9 (x 


-y) 


= 9 1 (x)9 2 (x)9 2 (y) 2 


- 9 (x) 2 9 (y)9 1 (y), 


(3.12) 


9 (0) 2 9 l (x + y)9 (x 


-y) 


= 9 (x)9 l (x)9 l (y) 2 


- 9 2 (x) 2 9 Q (y)9 2 (y), 


(3.13) 


9 (0) 2 9 2 (x + y)9 (x 


-y) 


= 9 (x)9 2 (x)9 Q (y) 2 


- 9 l (x) 2 9 1 (y)9 2 (y), 


(3.14) 


9 (0) 2 9 (x + y)9 l (x 


-y) 


= 9 (x)9 l (x)9 (y) 2 


-9 2 (x)%(y)9 2 (y), 


(3.15) 


9 (0)%(x + y)9 l (x 


-y) 


= 9 (x)9 2 (x)9 2 (y) 2 


- 9 1 (x) 2 9 (y)9 l (y), 


(3.16) 


9 (0) 2 9 2 (x + y)9 l (x 


-y) 


= 9 l (x)9 2 {x)9 { {y) 2 


- 9 (x) 2 9 (y)9 2 (y), 


(3.17) 


9 (0) 2 9 (x + y)9 2 (x 


-y) 


= 9 (x)9 2 (x)9 l (y) 2 


- 9 1 (x) 2 9 (y)9 2 (y), 


(3.18) 


9 (0) 2 9 l (x + y)9 2 (x 


-y) 


= 9 l (x)9 2 (x)9 (y) 2 


- 9 (x) 2 9 l (y)9 2 (y), 


(3.19) 


9 (0) 2 9 2 (x + y)9 2 (x 


-y) 


= 9 (x)9 1 (x)9 2 (y) 2 




(3.20) 



We give the proof of Lemma [331 in the appendix. 
Proof of Proposition 13.21 

The duplication formulas (13.111) are obtained by putting x = y = z in eqs. (13.121) , (13.131) and (13.141) . 
In order to prove eqs. (13.91 ) and (13.101 ), we first note that it follows by definition that 

9 (-z) = -0!(z), 9 2 (-z) = -9 2 (z), (3.21) 

and hence 

9 l (0) = -0„(O), 2 (O) = 0. (3.22) 



iA'(0,r) 

(3.10) 



(3.11) 
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From Lemma [331 we obtain the addition formulas for (p(z) and iff(z) (see eq. (13.101) ) as 



(p(x)(f(y) 2 - if/(x) 2 i//(y) 
(p(x)iff(x)iff(y) 2 - (f(y) 
iff(x)ifr(y) 2 - <p(x) 2 <p(y) 
(fix) - ifj{x) 2 ip(y)ifj(y) 



<p(x)if/(x) - <p(y)iff(y) 

<p(x + y) = , (3.25) 

i//(x)(p(yY - (p(xyif/(y) 



if/(x) - tp(x) 2 tp(y)if/(y) 
ip{x)ifj{x)ifj{y) 2 - (f(y) ' 
ip(x)ifj(x)(p(y) 2 - ijAy) 
(p(x) - i//(x) 2 <p(y)i//(yy 
ip{x)ifj(y) 2 - tff(x) 2 <p(y) 
ifj(x)y(y) 2 - <p(x) 2 i{f(y) ' 



*K*+y) = , ; (3 - 26) 

«*<.Y + v) = ^-w y >-.^) 2 

' v J 7 _ / _ _\ I/..A2 /-A.I./.A' 



<K* + ;y) = , — yv y . (3.28) 



Differentiating eqs. (|3.23l) and (13.251) by y and putting y = 0, we have 



<p (x) = -<p (0)ip(x) - ifr (O)iff(x) , <p(x) = — , (3.29) 



respectively. Here we have used 

cp(0) = -h <A(0) = 0, (3.30) 

which follows from eq. (13.221) . Equating the right hand sides of the two equations in eq. (13.291) , we 
have 

1 + (p(x) 3 + <AM 3 = -3^ (f(x)Hx), (3.31) 

iff '(0) 

which yields eq. (13.91 ) by multiplying #o(z) 3 . This completes the proof. □ 
Consider the map 

Cbz^ [6 (z) : 0i(z) : 6 2 (z)] e P 2 (C). (3.32) 

From the relations 

G k (z + 1) = -6 k (z), e k (z + r) = -e 3niT - 67Tiz e k (z) (k = 0,1,2), (3.33) 

we see that this induces a map from the complex torus L T = C/(Z + Zt) to Zy which is known 
to give an isomorphism L T (see, e.g. [2]). Since h [1 : -1 : 0], the addition formulas 

(I3.12D - (|3.20D induce the group structure on E T with the origin [1 : -1 : 0]. Denoting the addition 
of two points [x : X\ : x 2 ] and [x' : x[ : x' 2 ] as [x : x\ : x 2 ] © fx,', : x\ : ^], eqs. (I3.12l) - (|3.20l) 
imply 



[xq : Xi : x 2 ] © [Xq : x\ : JCj] 

r '2 2 i i i2 2ii i2 2ii-i 1A\ 

XqX\X\ — X 2 XnX 2 ■ XqX 2 Xq — X^X^X 2 \ \J.JH-) 



[XqX^Xq X 2 X^X 2 
\_XqX 2 X^ — X[XqX 2 



i2 _ 2 i i , i2 _ 2 i 1 1 /o ojr\ 

-^0-^2-^2 X|XqX| . X\X 2 Xy X^XqX 2 \ 



i2 2 i i . i2 2 i i n /-n 

X\X 2 Xq — XqX^X 2 . XqX\X 2 — X 2 XqX[\. {j.jO) 
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In particular, when the two points are equal, the duplication formula is given by 

2[xo : X\ : X2] = [xiC*5 — x o) • x o( x \ ~ x l) '■ x 2(^o ~ (3.37) 
Moreover, the inverse of [xo : x\ : ^2] is given by 

— [xq : X\ : Xq\ = \x\ : Xq : X2]. (3.38) 

We finally remark that jx can be also expressed as follows. Differentiating both equations in 
eq. (13.101) and putting z = 0, we have by using eq. (13.221) 

0'(O) 6'{0) 

«' (0) = 2 im- *' (0) = t(fly (3 ' 39) 

which yield 

^'(0) / (0) 

4 Ultradiscretization 

So far we have constructed the piecewise linear map (12.81) as the duplication map on the tropical 
cubic curve C K , whose general solution is given by eq. (12.321) . We have also presented the rational 
map (13.51) which arises as the duplication map on the Hesse cubic curve E^. The general solution 
of the map is given by eq. (13.81) . In this section, we establish a correspondence between the two 
maps and their general solutions by means of the ultradiscretization. 

4.1 Ultradiscretization of map 

The key of the ultradiscretization is the following formula: 

lim elogfe^ +e? +•-) = max[A,5, . . .]. (4.1) 

Putting 

x n = e~?, y„ = e~?, (4.2) 

we have from eq. (13.51 ) 

3r„+em \ 



(4.3) 



l+e « J + F„-elog(e « +e « J, 
l+e « J+X„-elog(e « +e * j, 

which yields, in the limit 6 — » +0, eq. (12.81) : 

Z n+1 = max[0, 3Z„] + F„ - max[3X„, 3Y n ], Y n+1 = max[0, 3F„] + X„ - max[3X„, 3F„]. (4.4) 

The limit of the invariant curve (|3.1I) yields C^: 

max[0,3X,3F] =X + F + K, (4.5) 
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by the use of 

3yu(r) = ei (4.6) 
In the above process of the ultradiscretization, we have calculated formally, for example, as 

elog(l -e 3J ^) = elog(l + e^) — > max[0,3X] (e -» +0). (4.7) 

However, when the original rational map contains the minus signs, such formal calculation some- 
times does not give consistent result. This may happen, for example, when we consider the limit 
of the exact solutions simultaneously, or when we consider the limit of the maps which are repre- 
sentation of certain group or algebra. In both cases, the cancellations caused by the minus signs 
play a crucial role on the level of rational maps, and the structure of the rational maps is lost be- 
cause such cancellations do not happen after taking the limit. This problem is sometimes called 
the minus-sign problem. 

Therefore, we usually consider the subtraction-free rational map to apply the ultradiscretization 
[45,46,49fl or we try to transform the map to be subtraction-free if possible [16]. Unfortunately, 
it seems that the map (13.51) cannot be transformed to be subtraction-free by simple transformations. 
However, in this case, it is possible to obtain valid ultradiscrete limit of the general solution in spite 
of the minus-sign problem. 

Remark 4.1 The nine inflection points of the Hesse cubic curve correspond to the vertices of the 
tropical cubic curve C K in the following manner; consider one of the inflection points [x : X\ : 
X2] = [1 : -1 : 0] = [gf : e^ 1 : e~?]. Then putting x t = ei (i = 0, 1,2) and taking the limit 
e — » +0, we have [X : X\ : X 2 ] = [0:0: -00] = [00 : 00 : 0]. Note here that on this level 
equivalence of the homogeneous coordinates is given by [X : X\ : X 2 ] = [X + L : X\ + L : X 2 + L] 
for any constant L. In the inhomogeneous coordinates, this point corresponds to (00, 00), which is 
linearly equivalent to the vertex V 3 = (K,K). Similarly, the two points [1 : -co : 0], [1 : -co 2 : 0] 
also correspond to V 3 . Furthermore, the triple of points {[1 : : — 1], [1:0: -to], [1:0: -co 2 ]} 
correspond to V 2 = (0, -K), and the triple {[0 : 1 : -1], [0:1: -to], [0:1: -co 2 ]} to Vi = (-K,0). 
In other words, three inflection points of the Hesse cubic curve degenerate to each vertex of Ck in 
the ultradiscrete limit. This explains the reason why the multiplicity of each vertex of C K is 3 and 
C K is not smooth while the Hesse cubic curve is non-singular. 

4.2 Ultradiscretization of general solution 

In this section, we consider the ultradiscrete limit of the solution. The following is the main result 
of this paper. 

Theorem 4.2 The general solution ({3.8\) of the rational map ({3.5\) reduces to the general solution 
d2.32D of the piecewise linear map < \2.8\) by taking the limit £ — > +0 under the parametrization 

X n = e^, F„ = A ~^ T = ~, zo = ^ll~\ u eR, K > 0. (4.8) 

5 It should be remarked that the term "tropical" has been used differently in the communities of geometry and 
integrable systems [18]. In the former community it has been used to mean piecewise linear objects, while in the 
latter subtraction-free rational maps. In the latter community the terms "crystal" or "ultradiscrete" have been used for 
piecewise linear objects. Therefore, it sometimes happens that the term "tropicalization" can be used with opposite 
meanings. 
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The ultradiscrete limit of the theta function can be realized by taking Im r — » 0, however, the 
limit of the real part of r should be carefully chosen in order to obtain consistent result [32]. For 
choosing the limit of real part of r, the following observation on the correspondence between the 
zeros of the theta functions and non-smooth points of S(u; a,fi, 6) is crucial. 

Observation: The ultradiscrete theta function 0(a; 6) defined by eq. (12.291 ) is a piecewise quadratic 
function with the period 1, and has zeros at u = n e Z. 0(w; 6) can be obtained from # (z;t) = 
#m T ) by taking the limit t — » [32]. Since the zeros of &o(z', r) are located at z = (m + pT + n 
(m, n 6 Z), the real zeros of &o(z; t) survive under the limit, giving the zeros of 0(w; 6) at u = n. 
From the definition of S given in eq. (12.301) and Figj6[ it is easy to see that the valleys at u = na and 
the peaks at u = fi + na (n e Z) of S correspond to the zeros of 0(^; 6) and @C—, 0\ respectively, 
as illustrated in FigQj)! In other words, valleys and peaks of the ultradiscrete elliptic function S 
arise from the zeros and poles of the corresponding elliptic function, respectively. Now, noticing 



i + na: zeiosof©(^ — 




Figure 10. Zigzag pattern of S(u; a, yS, 9) = 0(f ; 6) - 0(^; 6) and the zeros of 0(f ; 6), 0(^; 6). 

that the zeros of &( a ,b)(z, r) are located at z = (-a + m + j)t + (-b + n + j) (m, n 6 Z), the zeros of 
8 k (z, t) (k = 0, 1, 2) are given by 

(z,t): z = \m+^jT+^(n-l), 

0i(z,t): z = |m + ijr+i(n-l), 
1 



6 2 (z,t): z = rm + -(n - \), 



respectively. It is obvious that the zeros and poles of x n = and y n = cancel each other, 
respectively, in the limit r — » 0, which yields trivial result. Let us choose r — » -|. Then the zeros 
of 6i(z, r) (i = 0,1,2) become 

Z 2 Z 1 
Sofcr): z= --- = - + -, 

Z 1 Z 2 
ftfcr): z=--- = - + -, 

Z 

6» 2 (z,r): 2=3' 

respectively, which give the zigzag patterns in the limit as illustrated in FigHU These patterns 
would coincide with the those in Figj5] after an appropriate scaling. 
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Figure 1 1 . Zigzag patterns obtained by the limit r 



Before proceeding to the proof of the Theorem 14.21 we prepare the modular transformation of 
the theta function, which is useful in taking the limit of r. 



Proposition 4.3 [9,28] 

K- m (cr ■ z, o- t) = e ni((rz)cz (cr + d)s e 2 ^"^ # m (z, r), 



(4.9) 



where 



in = (///]. ///t). cr = | ^ ^ ) e SL:(Z), 



cr ■ m = mcr 1 + -(cd, a/?), cr • r = 



1 



ar + & 
cr + J' 



cr • z = 



CT + d 



(4.10) 



0m(cr) = _ ^ [^ m i + flcm 2 _ 2bcm l m 2 - ab{dm\ - cm 2 )\, 
k{o~) : an eighth root of 1 depending only on cr. 
Remark 4.4 Explicit expression of /c(cr) is given by [8] 



k(ct) 



Y flfa'fl 1 i acd^ c\ ( & 

" |c 



c 

I \\d\ 



c : odd 

^^(sgn^-lXsgnW-l)^ g f c . eyen 



where is the Legendre (Jacobi) symbol for the quadratic residue. In particular, (jj = 1. 



Proof of Theorem 14.21 The first key of the proof is to apply the modular transformation on 
9k(z, t) = |)(3z, 3r) (fc = 0, 1, 2) specified by 



cr 



1 

1 1 



From 



'k 7 3\ 9 
cr • ^ = (- - -,-), 0,«(cr) = --, k(<t) = 1 , 



(4.11) 



(4.12) 
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we have 



We put 



6 k (z,T) = 0, 4 i ^(3z,3t) = e (3t+ e* m 0,* 



T T = --, e>o, 



V T+ 3 T+ 3, 



T + 



(4.13) 



(4.14) 



and take the limit of e — > +0 which corresponds to r — > -|. Noticing = 3(^ + l)z, we have 

* (H>i) [ _L_] = £ e -fc-H> 2 +f M-l* ^HX**). (4 . 15) 

VT+ 3 T + 3 j ngZ 

The second key is to use the freedom of the imaginary part of z e C. Putting 



Z = — + iv, K, u, v 6 R, 
yA 



eq. (14.151) is rewritten as 



e f 



[*-»-!+?] 2 -6^+|-?)v ^i(n+|-i)[|v+^+|] 



f/tL 



If u and v satisfy 



8 u n 



then eq. (14.171) is simplified as 



&<k_7 3, 

( 1 6'2> 



-67r(n+4-J)v „3nin 



8\u-K(k+\) 



Therefore, asymptotic behaviors of x n and y„ as 6 ~ +0 are given by 



y n = 



6* 2 fc t) g 27Ti ^ e-?[ 2 #-5-("- 2 )] 2 e3n m - 

neZ 

8\{z,t) „ eZ 



(no+n2)jri 



e -§K«)H] : 



neZ 



-?[((^)H] 2 



e -f[((^)Hr 



(4.16) 



(4.17) 



(4.18) 



(4.19) 
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respectively, where 



n = Floor 



u-K 



3K 



+ 2, n x = Floor 



u-2K 
3K 



+ 2, n 2 = Floor 



u-3K 
3K 



+ 2. 



(4.20) 



Here, we note that z and u are understood as 2"zq and 2 n u Q , respectively. 

In order to obtain the final result, the remaining task is to relate the parameters 9 with K in the 
limit e — » +0. This can be done by considering the limit of the conserved quantity //(r) given by 
eq. (13.91 ) or eq. (13.401) . Noticing eq. (|3.40l) . we differentiate 8q(z, t) and 6 2 (z, r) by z after applying 
the modular transformation. Then 6' (0, r) and 6{(0, r) can be calculated by using eqs. (14.131) and 
(l4~T5b as 



6'(0,t) = h (3z,3r) 



dz 



z=0 



= ± e -<^y ( i + JLY e \* y e -\ («-d 2 +? ^s^j) 

dz \ in el ^ 

neZ 

-('♦=)'*SH* + 3H) 

«eZ L N ' 



z=0 



0'(O,r) = 3,(3z,3t) 



dz 



z=0 



= ± e -<i+&y ( i + JLY e & y e ->-^ + ^ n 

dz \ inel ^ 



\)z 6m(n-i)(z+i) 



neZ 

respectively, which imply 



.-=o 



3, = -6^=-6 



e' 2 (0,r) c i«*W„ _!)«-!(» 

neZ V ' 



;) 2 +3rari 



2 C T 2 e 
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as 6 — > +0. Accordingly, from eq. (|4.6I) we may put consistently as 



(4.21) 



Let us set x n = e " and y„ = e^" in eq. (14.191) . Then the complex factors in eq. (14.191 ) disappear in 
the limit of e — > and we finally obtain by using eq. (14.211) 



x - 9 -L 

2 



u - 2K\\ 1 
3£ ]/ 2 



+ 



9K 



u-3K\\ 1 



+ 



9K 



3K )) 2 
u - 3K\\ 1 
3K 1/2 



(4.22) 
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which is equivalent to eq. (|2.32l) . This completes the proof. 



□ 



We finally remark that the choice of parametrization (14.81) can be also justified by the following 
observation. The asymptotic formula (14.191 ) shows that (x n ,y n ) is in K 2 and that the quadrant of 
(x n ,y n ) changes according to the value of ((^)) as described in Tabled] Note that (x n ,y n ) never 
enters the first quadrant. It implies that qualitative behavior of the real orbit of the map (12.81) 
discussed in Section 3 is preserved under the limiting process. 



((#)) 


n n\ n 2 


n + n 2 n,\ + n 2 + 1 




[f,D 


N N N-\ 


2N - 1 2N 


(- +) 


L 3' 


N N-\ N - 1 


2N -I 2N -I 


(-, -) 


[0,0 


N- 1 N- 1 N- 1 


2N-2 2N-1 


(+,-) 



Table 1. Quadrant of (x n ,y n ) for e ~ +0, where N = Floor ((-^)) + 2. 
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A Proof of Lemma 3.3 



In this appendix, we give a proof of Lemma [331 Besides the cases k = 0, 1 , 2 we also use the cases 
k = |, |, | as well. Note that 6 k+3 (z, r) = 6 k (z, t). 

First, we remark that 6i(z,r) = ^ (0> 3)(3z, 3t) and 9 2 {z,t) = # ( i 3)(3z, 3t) are &(3z, 3r) and 
■&i(3z, 3r) of Jacobi's notation, respectively. Let us start from the eq. (A)-(4) in [14]: 

0i(w)0i(x)0i(y)6/i(z) - e 2 {w)e 2 {x)e 2 {y)6 2 {z) 

= ei(w')e k (x')9i(y')ei(z') - e 2 (w')e 2 (x')e 2 (y')e 2 (z'), iAA) 

where 

, w + x + y + z , w + x-y-z , w-x + y-z , w-x-y + z 

W = 2 ' X = 2 ' y = 2 ' Z = 2 • (A - 2) 

Replacing w as w — » w + |, we have w' — > W + |, x , —¥xf + ^,y f —*y' + ^ and z' — » z' + |. By 
using the formula 

0* (z + = g-^-fe g t+ i (z), 6 k (z + l) = e-™ e'^~ 2 ^ 6 k+i (z), (A.3) 
which easily follows by definition, we obtain 



- ei_(w)e h {x)e h (y)e h (z) + e Q (w)e 2 (x)e 2 ( y )e 2 (z) 
= ^(w')^ Wi0/)#i(z') - edw')e 5 _{x')6 { (y')e^z'). 
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(A.4) 



Replacing further as w — > w + |, x— >x+|,y— >y + | and z —* Z + f , we see that w' — > w' + y 
and x',y' , z' are unchanged. Then we obtain 

- 0|(w)03(x)0|(y)0|(z) + 0i(w)0 o (x>0 o (y)0o(z) 

= 9 (w')9 1 (x>)9 l ( y ')9 1 (z') - e^e^'w^yw^z'). iAS) 

Application of the same transformation to eq. (|A.1I) yields 

ei(w)6i(x)9i(y)6i(z) - 9 (w)9 (x)9 (y)9 (z) 

(A.6) 

= 9s_(w')9,{x')9 k {y')9 h {z') - 9 l (w')9 2 (x')9 2 (y')9 2 (z'). 

We put w = -(x + y + z) in eqs. (IA.4I ) and (IA.5I ). Then w' = 0, x' = -(y + z), y' = -(z + x) and 
z' = -(x + y). By definition it follows that 

9 (-z) = -0!(z), 0i(-z) = 0i(z), 2 (-z) = -0 2 (z), 0|(-z) = 0|(z), (A.7) 

and hence 

0j(O) = -0 O (O), 95(0) = 03(0), 2 (O) = O. (A.8) 
Therefore eqs. (IA.4I) and (IA.5I) yield 



- 05 (x + j + z)0i (x)0i (y)0i (z) -9 1 (x + y + z)9 2 (x)9 2 (y)9 2 (z) 

= 0i(O)0 o (y + z)0 o (z + x)0 o (x + y) - 0|(O)03 Cv + z)0|(z + x)03 (x + y), 

03 (x + y + z)0|(x)0| (y)03(z) + o (x + y + z)0 o (x)0 o (y)0 o (z) 
= -0 o (O)0 o (y + z)0 o (z + x)9 (x + y) - 03 (0)03 (y + z)03(z + x)03(x + y), 



(A.9) 



(A.10) 



respectively. Similarly, putting w = x+y+z, we have that w' = x+y+z and x' , y' , z' are unchanged. 
Then eq. (IA.6I ) yields 

03 (x + y + z)0|(x)03(y)03(z) - 9 (x + y + z)0 o (x>0 o (y)0o(z) 

(A.ll) 

= 9 s _(x + y + z)9i 2 (x)0i (y)0i (z) - 9 x {x + y + z)9 2 (x)9 2 (y)e 2 (z). 

Then from - (1A.9I )+ (1ATT01) + (1A.1 II) and dividing it by 2, we have 

o (O)0 o (y + z)0o(z + x)0o(x + y) = o (^ + y + z)0o(x)0o(y)0o(z)-0i(x+y + z)0 2 (x)0 2 (y)0 2 (z), (A. 12) 

which yields eq. (13.121) by putting z = -y in eq. (IA.12I) . Other addition formulas are derived from 
eq. (13.121) . Applying x — » x + |, y — » y + | on eq. (13.121) . we have 

o (O) 2 0!(x + y)0 o (x - y) = o (x)0 1 (^(y) 2 - 2 (x> 2 o (y)0 2 (y). (A.13) 
Repeating the same procedure on eq. (IA.13I) , we obtain 

o (O) 2 2 (x + y)9 (x - y) = 9 (x)9 2 (x)9 (y) 2 - 0i(x> 2 0i(y)0 2 (y). (A.14) 
Exchanging x <-» y in eq. (|3.12l) . we have 

0o(O) 2 o (x + y)0!(x - y) = o (x)0i (x)9 (y) 2 - 9 2 (x) 2 9 l (y)9 2 (y). (A.15) 
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(A.16) 
(A.17) 



respectively. Applying x —* x + 5 to eq. (IA. 171) . we obtain 

e (0) 2 e (x + y)e 2 (x -y) = 9 (x)e 2 (x)e 1 (y) 2 - e^xfe^e^). 



(A. 18) 




(A.19) 
(A.20) 



respectively. This completes the proof of Lemma [3731 □ 
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